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Abstract 

We demonstrate that the energy spectra of Ultra High Energy Cosmic 
Rays (UHECR) as observed by AGASA, Fly's Eye, HiRes and Yakutsk detec- 
tors, have the imprints of UHE proton interaction with the CMB radiation as 
the dip centered at E ~ 1 x 10 19 eV, beginning of the GZK cutoff, and very 
good agreement with calculated spectrum shape. This conclusion about proton 
composition agrees with recent HiRes data on elongation rate that support the 
proton composition at E > 1 x 10 18 eV. The visible bump in the spectrum at 
£~4x 10 19 eV is not caused by pile-up protons, but is an artifact of multiplying 
the spectrum by E 3 . We argue that these data, combined with small- angle clus- 
tering and correlation with AGN (BL Lacs), point to the AGN model of UHECR 
origin at energies E < 1 x 10 20 eV. The events at higher energies and the excess 
of the events at E > 1 x 10 20 eV , which is observed by AGASA (but absent in 
the HiRes data) must be explained by another component of UHECR, e.g. by 
UHECR from superheavy dark matter. 

1. Introduction 



The nature of signal carriers of UHECR is not yet established. The most 
natural primary particles are extragalactic protons. Due to interaction with the 
CMB radiation the UHE protons from extragalactic sources are predicted to have 
a sharp steepening of energy spectrum, so called GZK cutoff (Greisen 1966, Zat- 
sepin, Kuzmin 1966). For uniformly distributed sources, the GZK cutoff is charac- 
terized by energy E\/2, where the integral spectrum calculated with energy losses 
taken into account becomes twice lower than the power-law extrapolation from 
low energies (Berezinsky, Grigorieva 1988); Ex/ 2 = 5.7 x 10 19 eV. 
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Fig. 1. The HiRes data (Sokolsky 2002) on mass composition (preliminary). The 
measured x max at E > 1 x 10 18 eV (triangles) are in a good agreement with 
QGSJet-Corsika prediction for protons. 



There are two other signatures of extragalactic protons in the spectrum: 
dip and bump (Hill and Schramm 1985, Berezinsky and Grigorieva 1988, Yoshida 
and Teshima 1993, Stanev et al 2000). The dip is produced due p + 7cmb — > 
p + e + + e~ interaction at energy E ~ 1 X 10 19 eV. The bump is produced by 
pile-up protons which loose energy in the GZK cutoff. As it was demonstrated 
by Berezinsky and Grigorieva (1988) (see also Stanev et al 2000), the bump is 
clearly seen from a single source at large redshift z, but it practically disappears 
in the diffuse spectrum, because individual peaks are located at different energies. 
We shall demonstrate here that what is seen now in the observed spectrum as a 
broad bump is an artifact caused by multiplication of the spectrum to E 3 . 

As we shall demonstrate here, the observed spectra by AGASA, Fly's Eye, 
HiRes and Yakutsk arrays have the imprints of extragalactic protons in the form 
of the dip centered at E ~ 1 x 10 19 eV and of the beginning of the GZK cutoff, 
and in the form of good agreement between predicted and observed spectra. The 
measurement of the atmospheric height of EAS maximum, x max , in the HiRes 
experiment (for preliminary data see FigEJ) gives the strong evidence in favor of 
pure proton composition at E > 1 x 10 18 eV. Yakutsk data also favor the proton 
composition at E > 1 x 10 18 eV (Glushkov et al 2000). 
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At what energy the extragalactic component sets in? 

According to the KASCADE data (Kampert 2001,Hoerandel 2002 ), the 
spectrum of galactic protons has a steepening at E « 3 x 10 15 eV (the first knee), 
helium nuclei - at E m 6 x 10 15 eV, and carbon nuclei - at E « 1.5 X 10 16 eV. It 
confirms the rigidity-dependent confinement with critical rigidity R c = E c /Z m 
3 x 10 15 eV. Then galactic iron nuclei are expected to have the critical energy 
of confinement at E c ~ 1 x 10 17 eV, and extragalactic protons can naturally 
dominate at E > 1 x 10 18 eV. This energy is close to the energy of the second 
knee (Akeno - 6 x 10 17 eV, Fly's Eye - 4 x 10 17 eV, HiRes - 7 x 10 17 eV and 
Yakutsk - 8 x 10 17 eV). In Fig|T]the second knee is shown according to Fly's Eye 
observations. It illustrates a possible transition to extragalactic cosmic rays. 
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Fig. 2. The second knee from Fly's Eye data (Bird et al 1994). 

The AGASA data (Hayashida et al 1999) which shows an excess of events 
from regions of Galactic Center and Cygnus at E ~ 10 18 eV further confirm 
the picture outlined above. Indeed, unconfined galactic particles (most notably 
protons) propagate quasi-rectilinearly and show the direction to the sources, while 
the diffuse flux is dominated by extragalactic protons. 

The model of galactic cosmic rays developed by Biermann et al (2003) also 
predicts the second knee as the "end" of galactic cosmic rays (iron nuclei) due to 
rigidity bending in wind-shell around SN. The extragalactic component became 
the dominant one at energy E ~ 1 x 10 18 eV (see Fig.l in Biermann et al 2003). 

The good candidates for the sources of observed UHE protons are AGN. 
They can accelerate protons up to energy E max ~ 10 21 eV (Biermann and Streit- 
matter 1987, Ipp and Axford 1991, Rachen and Biermann 1993), they have power 
to provide the observed flux of UHE protons (Berezinsky, Gazizov and Grigorieva 
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2002) and finally there is direct correlation (Tinyakov and Tkachev 2001, 2003 
and references therein) between directions of arrival of UHE particles with ener- 
gies (4 — 8) x 10 19 eV and directions to BL Lacs, which comprise some particular 
class of AGN. 

Does it mean that UHECR puzzle has been already resolved in most con- 
servative way? 

In this model AGN cannot be the sources of observed particles with energy 
E > 1 x 10 20 eV (Berezinsky, Gazizov, Grigorieva 2002): the attenuation length 
for a proton of this energy is smaller than 135 Mpc, and correlation with AGN 
would be seen for all these particles, contrary to observations. The particles 
observed at E > 1 x 10 20 eV, in particular those detected in AGASA, imply the 
presence of another component, e.g. produced by decays of superheavy DM. 

2. Calculation of the spectra: dip and GZK cutoff 

We calculate the UHE proton spectra in the model with uniform distri- 
bution of the sources, with CR luminosity of a source L p (z), and CR emissivity 
C(z) = n(z)L p (z) = C (l + z) m , where n(z) is comoving space density of the 
sources at epoch with redshift z, index refers to z — 0, and m describes cosmo- 
logical evolution of the sources. We assume the generation spectrum of a source 

Q g (E g , z) = J^ {Z J^ q gcn (E g ), (1) 

„ ( F )-l l / E l at E 9< E c (9) 

W ' J I E^(E g /E c )-^ at E g > E c [Z) 
The diffuse spectrum can be calculated as 

J P (E) = ^ ln ^° + ^ J Q ZmaX dt(l + zTq gcn (E g (E, z), E) ^L, (3) 

where E g (E,z) is generation energy of a proton at epoch z, dE g /dE is given in 
Refs. (Berezinsky, Grigorieva 1988, and Berezinsky, Gazizov, Grigorieva 2002), 
and dt is given as 

dz , lS 

dt = (4) 

h (i + z) y/n m (i + z y + n A 

with H Q , fl m , and fl\ being the Hubble constant, relative cosmological density 
of matter and vacuum energy, respectively. In FigET] we present the calculated 
spectra compared with AGASA (Takeda et al 2002), HiRes (Abu-Zayyad et al 
2002a), Fly's Eye (Bird et al 1994) and Yakutsk (Glushkov and Pravdin 2001) 
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Fig. 3. Comparison of calculated spectra for non-evolutionary model with experi- 
mental data. The normalization has been first performed by the AGASA data at 
E = 1 x 10 18 eV, adjusting the emissivity Cq. To fit the data of HiRes, Fly's Eye and 
Yakutsk the emissivity Cq has been scaled by factors 0.63, 0.8 and 1.7, respectively. 



data. The dip is well confirmed by all data. HiRes and Yakutsk data agree well 
with existence of GZK cutoff, while the AGASA data show the excess of events 
at E > 1 x 10 20 eV. 

For the curves above we used the following parameters: 7 g = 2.7, E c = 
1 x 10 18 eV, £ max = 1 x 10 21 eV, H = 70 km/s Mpc, Vt m = 0.3, fi A = 0.7. 
We have chosen the CR emissivity Cq = 3.5 x 10 46 erg/Mpc 3 yr to normalize the 
AGASA data, and scaled Cq by factors 0.63, 0.80 and 1.7 to normalize the data 
of HiRes, Fly's Eye and Yakutsk, respectively. 



3. Ein as characteristic of the GZK cutoff 

Ex/2 is the energy where the flux calculated with energy losses becomes 
twice less than power-law extrapolation of integral spectrum. In FigR~k the func- 
tion E^" 1 ^ J(> E) is plotted as function of energy (7 > 7 9 is the effective index). 
For wide range of generation indices 2.1 < 7 9 < 2.7 the cutoff energy is the same, 
E1/2 ~ 5.7 x 10 19 eV. We have determined £"1/2 from the Yakutsk data. For this 
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Fig. 4. E1/2 as numerical characteristic of the GZK cutoff. In panel a) the calculations 
for different 7 gen are presented. In panel b) Ei/ 2 is found from the integral spectrum 
of the Yakutsk array using two fits of the integral spectrum. 



we found two fits of the Yakutsk integral spectrum with help of trial functions, 
as shown in Figl4~b. They have good x 2 / n equal to 0.65 and 0.52. The corre- 
sponding values of E1/2, 5.6 x 10 19 eV and 6.2 x 10 19 eV, agree well with the 
theoretical value. Note, that in the fits above x 2 /n are the formal values from 
which probabilities cannot be calculated in the standard way, because the points 
in the integral spectrum are correlated quantities. 

This analysis obviously cannot be extended to the AGASA integral spec- 
trum, because of too many events at the highest energies. Unfortunately, we do 
not have the HiRes integral spectrum to perform the analysis as that above. 

4. Bump in the diffuse spectrum 

The bump is distinctly seen in the measured spectra when they are multi- 
plied to E 3 (as example see the HiRes spectrum in Fig lS~]) . 
Is this bump really composed by pile-up protons? 

To discuss the bump produced by pile- up protons it is convenient to intro- 
duce the modification factor rj (Berezinsky and Grigorieva 1988), defined for the 
power-law generation spectrum according to Eqs.fpQ) - (jBJ). 
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Fig. 5. HiRes spectrum from Abu-Zayyad et al 2002b. 



Without loss of generality, one can assume the power-law generation spec- 
trum, and present the diffuse spectrum J P (E) as the product of unmodified flux, 
J unm , and modification factor T)(E, z max ), which describes the distortion of spec- 
trum by energy losses: 



JuUE) = ^( lg -2)^-E-^, (5) 



— \-r°(E,z) dE g 
ri(E,z max ) = dz ^=—±, ( 6 ) 

(i + z)J(i + zfn m + n A dE 



and 



where X(E, z) = E g (E, z)/E. When E is small, A = (1 + z), dE g /dE = (1 + z) 
and modidification factor tends to energy-independent value 

r/(^ max ) = /_ dz - [ + ) _■ (7) 

[i + zfn m + n A 
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Fig. 6. Modification factors for diffuse spectra with different z max . The curves be- 
tween z max = 2.0 and z max = 0.2 have z max = 0.3 ,0.5 and 1.0. The solid curves 
are for 7 5 = 2.0 and dashed ones - for 7 9 = 2.7. The pile-up peaks, clearly seen at 
small z max , disappear when summation of sources goes to large 2 max > 1. 

In Fig ED the evolution of modification factor with growth of z max is shown. 
The pile-up peaks are seen at small z max . For large z max they are located at differ- 
ent energies, and their sum is given by a smooth curve without visible bump(s). 
The solid curves correspond to 7 g = 2.0 and the dashed ones - to j g — 2.7. In 
FigEJwe present two modification factors: r\ ee (dotted curve), when the adiabatic 
energy losses and that due to pair-production are included, and rjtot (solid and 
dashed curves), when all energy losses are included. Both modification factors are 
normalized at low energies by values r)(z max ), given by Eq.(JZJ), so that their low- 
energy limit is equal 1. In case of r] ee (dotted curve), it tends to 1 for both high 
and low energies, because pair-production energy losses disappear there. This 
curve is shown for the case 7 g = 2.7. 

Transition to the photopion energy losses is accompanied by practically 
invisible bump. 

If one displays the spectrum without factor E 3 , it also shows no presence 
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of the bump. 
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Fig. 7. Modification factor as characteristic of the dip and bump. The dotted curve 
show r] ee , when adiabatic and pair-production energy losses are included, for the 
case j g = 2.7. The solid and dashed curves include also the pion-production energy 
losses. The pile-up peaks are practically absent. 



5. AGN as UHECR sources 

AGN are traditional candidates for UHECR sources. The particles can 
be accelerated there up to -E max ~ 10 21 eV (Biermann and Streitmatter 1987, 
Ipp and Axford 1991, Rachen and Biermann 1993). The CR emissivity £ ~ 
3 x 10 46 erg/Mpc 3 yr is well within total emissivity of AGN, e.g. that of Seyfert 
galaxies is of order £$ y ~ risy-^sy ~ 1 x 10 48 erg/Mpc 3 yr. AGN origin of UHECR 
results in presence of the GZK cutoff, and this prediction agrees with data of 
Fly's Eye, HiRes and Yakutsk detector. In case AGASA spectrum is the correct 
one, the extra component at E > 1 X 10 20 eV is needed. In FigJHDwe plot in 
the AGASA panel the spectrum of UHE photons from Superheavy Dark Matter 
(SHDM), which is according to recent calculations (Berezinsky and Kachelriess 
2000, Sarkar and Toldra 2001, Barbot and Drees 2002) is Ij(E) oc E~ 2 . Note 
that in this case at energies E < 1 X 10 20 eV the protons dominate. 

The remarkable direct evidence for AGN as UHECR sources at energy 
range (2 — 8) x 10 19 eV has been found in series of papers by Tinyakov and 
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Tkachev 2001,2003 and references therein. There is the correlation between arrival 
directions of UHE particles in the AGASA and Yakutsk detectors and directions 
to BL Lacs, which are AGN with jets directed towards us. This correlation implies 
that UHE signal carriers propagate rectilinearly in the universe. The rectilinear 
propagation of the signal carriers from the point-like sources gives also the most 
natural interpretation of the small-angle clustering (Dubovsky, Tinyakov, Tkachev 
2000; Fodor and Katz 2000) observed by AGASA (Uchihori et al 2000). 

If UHECR are protons, as we argue above, the extragalactic magnetic 
fields should be very weak. 

6. Extragalactic magnetic fields and rectilinear propagation of UHE 
protons 

How weak the magnetic field must be ? 

Magnetic field must not produce the angular deflection larger than an- 
gular resolution of sources in the detectors, which is typically 9 TCS « 2.5°. The 
correlation is found in the energy range (4 — 8) x 10 19 eV, for which the largest 
attenuation length is Z att ~ 1000 Mpc. The required upper limit for the magnetic 
field, which is homogeneous on this scale, is Bi < 2 x 10~ 12 1^q 0Q G, where /i oo is 
attenuation length for 4 x 10 19 eV protons in units of 1000 Mpc. For a magnetic 
field with small homogeneity length, Itom, the required upper limit is 

B < E ° rc * ~ 6 x 10- 10 G, (8) 

'att'hom 

where the numerical value is given for l att ~ 1000 Mpc and Zhom ~ 10 kpc. 
We argue that these fields are not excluded. 

The observed Faraday rotations give only the upper limits on large scale 
extragalactic magnetic field (Kronberg 1994). All known mechanisms of gen- 
eration of the large scale cosmological magnetic field results in extremely weak 
magnetic field ~ 10~ 17 G or less (for a review see Grasso, Rubinstein 2001). The 
strong magnetic field can be generated in compact sources, presumably by dy- 
namo mechanism, and spread away by the flow of the gas. These objects thus 
are surrounded by magnetic halos, where magnetic field can be estimated or mea- 
sured. The strong magnetic fields of order of 1/xG are indeed observed in galaxies 
and their halos, in clusters of galaxies and in radiolobes of radiogalaxies. As an 
example one can consider our local surroundings. Milky Way belongs to the Local 
Group (LG) entering the Local Supercluster (LS). LG with a size ~ IMpc contains 
40 dwarf galaxies, two giant spirals (M31 and Milky Way) and two intermediate 
size galaxies. The galactic winds cannot provide the appreciable magnetic field 
inside this structure. LS with a size of 10 - 30 Mpc is a young system where 
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dynamo mechanism cannot amplify significantly the primordial magnetic field. 
In fact LS is filled by galactic clouds submerged in the voids. The vast ma- 
jority of the luminous galaxies reside in a small number of clouds: 98 % of all 
galaxies in 11 clouds (Tully 1982). Thus, accepting the hypothesis of generation 
of magnetic fields in compact sources, one arrives at the perforated picture of 
the universe, with strong magnetic fields in the compact objects and their halos 
(magnetic bubbles produced by galactic winds) and with extremely weak mag- 
netic fields outside. However, even in this picture there is a scattering of UHE 
protons off the magnetic bubbles and the scattering length is l sc ~ l/7rR 2 n, where 
R is the radius of magnetic bubbles and n is their space density. Among various 
structures, the largest contribution is given by galaxy clusters which can provide 
l sc ~ (1 - 2) x 10 3 Mpc. 

7. Discussion 

In this paper we consider AGN as UHECR sources. We assume that there 
are two mechanisms of acceleration operating there: one accelerates protons along 
the jet (Chen et al 2002), while the other operates in the shock where the jet 
terminates (Biermann and Streitmatter 1987, Ipp and Axford 1991, Rachen and 
Biermann 1993). When a jet is directed towards us, we observe the jet accelerated 
protons and classify the AGN as BL Lac. The shock-accelerated protons are 
emitted isotropically, and their flux at the earth is smaller than that from the jet 
component because of the geometrical factor. We assume additionally that CR 
luminosity of the jet component is larger than of isotropic one. Then BL Lacs 
are to be observed as UHECR sources at large distances. At small distances, 
less than 135 Mpc (attenuation length of proton with energy 1 x 10 20 eV), the 
total number of AGN is smaller, and the probability to observe an AGN as BL 
Lac is correspondingly smaller. But AGN can be detected at these distances due 
to isotropic component. At the distances less 31.5 Mpc (attenuation length of a 
proton with energy 2 x 10 20 eV) there are as minimum 11 Seyfert galaxies (with 
redshifts z < 0.009) and 4 nearby radiogalaxies. The lower limit on the number of 
UHECR sources for this distance is estimates by Tinyakov and Tkachev 2003 on 
basis of small-angle clustering as 30 at 90% CL and 4 at 99% CL, in agreement 
with observed number of AGN given above. For 135 Mpc (E p = 1 x 10 20 eV 
the number of sources become 80 times larger. It justifies the use of uniform 
distribution of the sources for flux calculations at E < 1 X 10 20 eV. 

For events at energies E > 1 x 10 20 eV the AGN model meets the difficulties 
due to short attenuation length of protons (31.5 Mpc for E p = 2 x 10 20 eV). 
UHECR sources must be seen in the direction of each proton with such energy. 
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One can consider two cases. 

(i) The AGASA data at E > 1 x 10 20 eV are included in analysis. The AGASA 
excess should be explained by another component as shown in Fig]3~l 

(ii) If to neglect AGASA data, one has to analyse only three events with energies 
E > 1 x 10 20 eV, namely one Fly's Eye event with Ek 3 x 10 20 eV, one HiRes 
event with E w 1.8 x 10 20 eV, and one Yakutsk array with E w 1.0 x 10 20 eV. 
Let us consider the highest energy Fly's Eye event. If it is induced by a proton 
and its source is an AGN, then the distance to it should be only 20 - 30 Mpc. No 
AGN is observed in this direction and at this distance. As in case (i) we again 
have to assume the new component at E > 1 x 10 20 eV, but this time it is based 
on very low statistics. 

To avoid the problems with E > 1 x 10 20 eV within the AGN model, one 
can assume the presence of strong extragalactic magnetic field. The AGN model 
in this case explains all data, except the AGASA excess and correlation with 
BL Lacs. In particular the small-angle clustering can be explained by magnetic 
lensing (Harari et al 1999, Sigl et al 1999, Yoshiguchi et al 2002). One can try 
to explain the correlation with BL Lacs by flux of light neutral hadrons from 
AGN (Berezinsky et al 2002, Kachelriess et al 2003). However, in the case we 
want to explain the dip observed in UHECR by primary protons, the flux of 
light hadrons produced in AGN must be subdominant with more flat spectrum 
(oc E~ 2,7 ) than protons. As a result the correlation must increase with energy, 
which is not observed. The number of correlated events is expected to be small, 
being proportional to the fraction of events induced by light neutral hadrons. 

8. Conclusions 

The preliminary data of HiRes indicate to proton composition of UHECR 
at E > 1 x 10 18 eV (Sokolsky 2002). The observed energy spectra reveal the 
signatures of interaction of UHE protons with the CMB in the form of the dip, 
beginning of the GZK cutoff and the good agreement with the predicted spectrum. 
Combined with small-angle clustering (Uchihori et al 2000) and correlation with 
BL Lacs (Tinyakov and Tkachev 2001, 2003), these data require the rectilinear 
propagation of UHE protons. The correlation with AGN (BL Lacs) becomes thus 
the most sensitive method of measuring extragalactic magnetic fields, in case they 
are very weak. 

Events with E > 1 x 10 20 eV imply the new component, which can be 
e.g. gamma-rays from the decay of superheavy particles. This case is especially 
favorable for 12 AGASA events with theses energies. The spectrum of UHE 
photons from the decay of superheavy particles fits well the AGASA observations 
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(see FigEO. 
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